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Abstract 

Vibration acceleration levels on large space platforms exceed the requirements of many 
space experiments. The Microgravity Vibration Isolation Mount (MIM) was built by the Canadian 
Space Agency to attenuate these disturbances to acceptable levels, and has been operational on the 
Russian Space Station Mir since May 1996. It has demonstrated good isolation performance and 
has supported several materials science experiments. The MIM uses Lorentz (voice-coil) magnetic 
actuators to levitate and isolate payloads at the individual experiment/sub-experiment (versus rack) 
level. Payload acceleration, relative position, and relative orientation (Euler-parameter) 
measurements are fed to a state-space controller. The controller, in turn, determines the actuator 
currents needed for effective experiment isolation. This paper presents the development of an 
algebraic, state-space model of the MIM, in a form suitable for optimal controller design. 
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Introduction 

Acceleration measurements on the U S. Space Shuttle and the Russian Mir Space Station 
show acceleration environments that are noisier than expected [1], The acceleration environment 
on the International Space Station (ISS) likewise will not be as clean as originally anticipated; the 
ISS is unlikely to meet its microgravity requirements without the use of isolation systems [1], [2], 
While the quasi-static acceleration levels due to such factors as atmospheric drag, gravity gradient, 
and spacecraft rotations are on the order of several micro-g, the vibration levels above 0.01 Hz are 
likely to exceed 300 micro-g rms, with peaks typically reaching milli-g levels [3], These 
acceleration levels are sufficient to cause significant disturbances to many experiments that have 
fluid or vapor phases, including a large class of materials science experiments [4], 

The Microgravity Vibration Isolation Mount (MIM) is designed to isolate experiments from 
the high frequency (>0.01 Hz) vibrations on the Space Shuttle, Mir, and ISS, while passing the 
quasi-static (<0.01 Hz) accelerations to the experiment [5], It can provide up to 60 dB of 
acceleration attenuation to experiments of practically unlimited mass [6], The acceleration- 
attenuation capability of the MIM is limited primarily by two factors: (1) the character of the 
umbilical required between the MIM base (stator) and the MIM experiment platform (flotor), and 
(2) the allowed stator-to-flotor rattlespace. A primary goal in MIM design was to isolate at the 
individual experiment, rather than entire rack, level; ideally the MIM isolates only the sensitive 
elements of an experiment. This typically results in a stator-to-flotor umbilical that can be greatly 
reduced in size and in the services it must provide. In the current implementation, the umbilical 
provides experiments with power, and data-acquisition and control services. Even with the 
approximately 70-wire umbilical the MIM has demonstrated good isolation performance [6]. 

The first MIM unit was launched in the Priroda laboratory module which docked with Mir 
in April 1996. The system has been operational on Mir since May 1996 and has supported several 
materials science experiments [1], An upgraded system (MIM-II) will be flown on the U S. Space 
Shuttle on mission STS-85 in July 1997 [6], 



In order to develop controllers for the MIM it is necessary to have an appropriate dynamic 
model of the system. The present paper presents an algebraic, state-space model of the MIM, in a 
form appropriate for optimal controller design. 

Problem Statement 

The dynamic modeling and microgravity vibration isolation of a tethered, one-dimensional 
experiment platform was studied extensively by Hampton [7], It was found that optimal control 
techniques could be effectively employed using a state-space system model, with relative-position, 
relative-velocity, and acceleration states. The experiment platform was assumed to be subject to 
Lorentz (voice-coil) electromagnetic actuation, and to indirect (umbilical-induced) and direct 
translational disturbances. 

The task of the research presented below was to develop a corresponding state-space model 
of the MIM. Translational and rotational relative-position, relative-velocity, and acceleration states 
were to be included, with the rotational states employing Euler parameters and their derivatives. 
The MIM dynamic model must incorporate indirect and direct translational and rotational 
disturbances. 


System Model 

A schematic of the MIM is depicted in Figure 1 . The stator, defined in reference 
frame is rigidly mounted to the orbiter. The flotor, frame ©, is magnetically levitated above 

the stator by eight Lorentz actuators (two shown), each consisting of a flat racetrack-shaped 
electrical coil positioned between a set of Nd-Fe-Bo supermagnets. The coils and the supermagnets 
are fixed to the stator and flotor, respectively. Control currents passing through the coils interact 
with their respective supermagnet flux fields to produce control forces used for flotor isolation and 
disturbance attenuation [5], 

The flotor has mass center F* and a dextral coordinate system with unit vectors f , f , and 
/ , and origin F 0 . The stator (actually, stator-plus-orbiter) has mass center S' and a dextral 

—3 

coordinate system with unit vectors | [ , $ 2 , and , and origin S 0 . The inertial reference frame © is 
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similarly defined by w,, « 2 , and and origin N 0 . The umbilical is attached to the stator at S u , and 
to the flotor at F u . 



Figure 1. Schematic of the MIM 


State Equations of Motion 

Translational Equations of Motion 

Let E be some flotor-fixed point of interest for which the acceleration is to be determined. 
If £ has inertial position r _ , then its inertial velocity and acceleration are f = — Ir I and 

~ N ° E “Vo E dt Wo Ej 

N d\ N d( \1 

Ln 0 e - —jf\ ~dt\^ E ) » res P ectlv ely- (The presuperscript indicates the reference frame of the 

differentiations. The subscripts indicate the vector origin and terminus.) The angular velocity and 
angular acceleration of the flotor with respect to the inertial frame are represented by N co F and 

N a F , respectively, where N a F = N o/) . 

Let F be the resultant of all external forces acting on the flotor; M_ Fi ' F (or simply A£), the 
moment resultant of these forces about m, the flotor mass; and lJ lF \ or / ), the central inertia 
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dyadic of the flotor for / , / , and / . Then Newton’s Second Law for the flotor can be 

—1 —2 —3 


expressed in the following two forms: 


E = ™t NoF - 


From Equation (2), 


M=l v+ 




It will be useful to find an expression for r Ni)E in terms of the acceleration r N ° s of the umbilical 
attachment point S u , and in terms of the extension of the umbilical from its relaxed position. 


Begin with the following: r 5A = r NoE + r EF . - r NoSo - r ^ . 


Differentiation of Equation (4) yields 


ts A = t NoE +V x LeFu - r NoS - Vx r SA 


(Eq. 5) 


A second differentiation gives 


tsA=^ 0 E+ V *L EF + % F x( Vxr ff J-r VA - V X^-Vxj VxrJ. (Eq. 6) 


Substitution for ^a^from Equation (3) into Equation (6) yields 


- V« ~ + 


L No e + |f"‘ • Mr N Q} f x (/• "a/)]} xr^-r 

- V s x ( V x L SoS ) + V x ( V x r^) . 


- V xr 

NqSq “ * Ls^ a 


(Eq. 7) 


In these equations 


0) = co + 0) . 


(Eq. 8) 


Under the assumptions that 0) s and N a s are negligibly small and , therefore, that 


-W 0 s 0 * > 


(Eq. 9) 



6 

Equation (7) reduces to 

£ s a ' Mr Q± x (| - ^ )JJ xf ef u ~ ^n 0 s, + S ^. Fx ( 5 ^. Fx ?!£f i ,)- (Eq. 10) 

Linearization about s g) F = 0 yields the following result: 

£ V, = + \l ' — } x —ef u ~ £v. - (Eq. 1 1 ) 

Appropriate expressions for F and M will now be determined, for substitution into Equations (1) 
and (11), respectively. Those amplified equations will then be used to obtain a more useful 
expression for r s F . [See Equations (43-48).] 

The force resultant F is the vector sum of the eight actuator (coil) forces F (/ = 1,...,8), 

with resultant of the umbilical force Fu,, caused by umbilical extensions from the relaxed 
position; of the direct disturbance forces, with resultant Fj, and of the gravitational force Fg. 

Gravity may be neglected for a space vehicle in free-fall orbit. The moment resultant M is the 
vector sum of the moments due to the coil forces, with resultant M c ; of the moment M_ ut due to 

the umbilical force F ut ; of the moment M ur due to the umbilical rotations from the relaxed 
orientation; and of the moment M_ d due to the direct disturbance forces. There is no moment due 

to gravity, since M is about the flotor center of mass F . In equation form, assuming the i' h coil 
force to be applied at point B h 

+£„+£* (Eq. 12) 

i=l 

8 

M = Tj r - F % x Ec + l rFm x F ut + M ur + M d 


and 


(Eq. 13) 
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More explicit expressions for F' c and F ut will now be developed. If the actuator has coil 

A A. 

current length L t , and magnetic flux density B i B i , then the associated actuator force 
becomes 

£ = /,A4x|, (Eq. 14) 


Assume a translational stiffness K' t for an umbilical elongation in the s direction, and a 

corresponding translational damping C,' . If a post-superscript u is used to refer to the umbilical in 
its relaxed (unextended, untwisted) position, then the force the umbilical exerts on the flotor 
becomes 



(Eq. 15) 


Define the following, for / = 1, 2, 3 : 

*oi = (-v^ ~ - vl) (Eq 16 ) 

and x bl =x ai . (Eq 17) 

If N qi s « 0 , Equation (15) becomes 


£„=-! 




(Eq. 18) 


L i=l J 

The relative positions x m and the relative velocities x bi will be six of the nine translational states used 
in the state-space formulation of the system equations of motion. 

As with F' c and F ul above, M ur can be also be expressed in more explicit form, in analogous 
fashion. Assume a rotational stiffness K' r and a rotational damping C ' r , for umbilical twist about the 
5 direction. Let <J ) F ^ s n, represent the rotation of the flotor, relative to the stator, from the relative 
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A 

position in which the / and s coordinate systems are aligned. Wfi is the rotation axis, and <{> is the 

angle of twist about that axis. Using the post-superscript u as before, the moment M ur can be 
expressed by the following: 


V F %) s]s i + ic r N J t {V' S n r V F, *$)-h 




(Eq. 19) 


or 


M.r = *‘ F/S «;) s,]|, + |:c[(* f ' s 2 < - f F %)-1 ]|,J. (Eq. 20) 

Equation (20) can be expressed in alternate form using Euler parameters. Let Fls h be 


described in © by 


' 7s /r=e,|, +e 2 s 2 +e 3 | 3 . (Eq. 21) 


Define the following Euler parameters [8], 


F/s P„=cos|, 

(Eq. 22) 

<^p, =^sin|. 

(Eq. 23) 

™%=e 2 »n|. 

(Eq. 24) 

< s « s P,=e,sin|, 

(Eq. 25) 

and F,i p=sin | F ' s h • 

1 

(Eq. 26) 

For small values of <)>, the Euler parameters can be simplified: 


F,s Mi, 

(Eq. 27) 

< w *P,=e,*/2, 

(Eq. 28) 
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<WS P ,=«,*/2, 

(Eq. 29) 



(Eq. 30) 

and 

=| F/S n ■ 

— L + 

(Eq. 31) 


Note that, for small angles, 

b F '%=2 F % (Eq. 32) 


This equation can be used to simplify the stiffness terms of Equation (20). 
For the damping term. Equation (3 1) can be differentiated to yield 



f/5 {3 = — cos — F/S n . + sin ^ F ! s ~ 
-2 2 “* 2 -4> 

(Eq. 33) 

or, for small angles, 

2 = <j) F,s rr + 4» F,s ir ■ 

— — 9 — 9 

(Eq. 34) 

Equation (20) now becomes 


M ur = -2 

t K‘, [( ■ ) ■ 1 ] l + 1 c;( • 1 ■) 1 , 

_ i=l 1 J i=l 

(Eq. 35) 

Define the following, for / 

= 1,2,3: 



x„ 

(Eq. 36) 

and 


(Eq. 37) 


The assumption that N co s is negligible yields, finally, 

+ (Eq. 38) 

Note that Equation (11) describes r s F in terms of the acceleration of an arbitrary flotor- 

fixed point E. For E located at flotor mass center F\ Equation (1 1) can be used straightforwardly 
with Equation (1) to yield 



10 


= “£ + (r (Eq. 39) 

Define now three unknown-acceleration terms, to be used with Equation (39). The first term 
represents the indirect translational acceleration disturbance input to the flotor, applied at the stator 


end of the umbilical: 


-in ~ ■ 


(Eq. 40) 


The second term represents the direct translational acceleration disturbance to the flotor, due to 


unknown disturbance force F_ d 


4 a =-F d . 

m 


(Eq. 41) 


And the third represents the direct angular acceleration disturbance input to the flotor, due to F, 


a^= /' ■ K d , 


(Eq. 42) 


Substitution from Equations (12), (13), (14), (18), (38), (40), (41), and (42) into (39) yields the 
following result: 

Zv . = +£„)+/■' (M, + M„ F +a d xr_ -a la +a d , (Eq. 4: 


£=!£=! Wi,xs,, 


1=1 i= 


£., = -Z(0.+c;xJ|, , 

- 1=1 


(Eq. 44) 


(Eq. 45) 


m c = Zm; = x (/mz, x I,) , 

/=! i=i 

3 

M u , = -r F . Fu x £(tfX+CX 1 )| j , 

- ,=I 


(Eq. 46) 


(Eq. 47) 
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and 


M„=-2|t[/:x+c>.)l. 

r=l 


(Eq. 48) 


Substitution from Equation (43) into Equation (11) produces the following equation for the 
acceleration of arbitrary flotor point E: 

t NaE = “ fcc + £«») + r • (Me + M ut Lf-e + x r rE + q d . (Eq. 49) 


Assuming % to be negligible, one also has the following: 



s d / \ 




(Eq. 50) 

and 



-V, “ V-vJ' 

(Eq. 51) 


(Note that assuming % s to be negligible does not imply that © and © are identical; it means 
rather that © can be treated as if it is in pure translation relative to © for the frequencies of 

interest.) Equations (43), (49), (50), and (51) provide the basis for a state-space form of the 
translational equations of motion, using x ai , x bi , and low-pass-filtered approximations to the s, 
components of r Nt>E [see Equations (96) and (103)], as states. 


Rotational Equations of Motion 

Let F/,s jp=sin ~ F/S n as before [Eq. (26)]. Differentiating the left side twice produces 
2 




(Eq. 52) 

and 

^ ( F/S p)-^( F/ 'p)+V x^( F,s v)+% s x F ' s p +V x'% 

(Eq. 53) 


Assuming as before that N co s ~ 0 , Equations (52) and (53) become, respectively, 
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(Eq. 54) 

and 

"*t=£m 

(Eq. 55) 

Similarly, 

F/5 p- F/s p“=^-( F /^-^“) 

(Eq. 56) 

and 

W |-"T =%r( Fl t- Fl t)- 

(Eq. 57) 


N d 2 


dt 


Returning to Equation (26), two differentiations of the right side yield 

( F/Sn \ ^ 4> F/S~ f<j>l • ^ FIS- , T d) . <t> F/S” 

\ §/ = 2 C0S I -fvfJ 2; cos 2 ^ +sin ? »#• 


(Eq. 58) 


Linearizing about <{) = 0 and <j) = 0 , and assuming F/S n^ to vary negligibly with time. Equation (58) 
becomes 2 F/ *p = <j> F/S n* . (Eq. 59) 

Similarly, Equation (33) reduces to 2 P = <j> F,s n^ . (Eq. 60) 

Equations (56), (57), (59), and (60) provide the basis for a state-space form of the rotational 
equations of motion, using as states the s, components of F/S p - f/s p" and of F/5 p - F/5 p“ (i.e., 
x a and x ej , respectively, for / = 1, 2, 3). 


Equations of Motion in State-Space Form 
From Equation (16), 

-V. _ = I +^a2^2 + *a3^3- (Eq. 61) 

Using Equations (16), (17), and (50), 

-V. — -V. = ^1 + X b2$-2 + 3 • 


(Eq. 62) 



Similarly, employing Equations (16), (17), and (51), 
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-v« — — V« = 1 + *i> 2^2 + *63^3 • (Eq. 63) 

Introduce the use of a presuperscript in parentheses to indicate the coordinate system used for 
componentiation. (This notation allows vectors to be expressed unambiguously in terms of their 
measure numbers.) Then Equations (61) and (62) take the respective forms. 



(Eq. 64) 



(Eq. 65) 


where x a and x b are defined as indicated. (Cf. Equations (16) and (17).) x d and x e have 
corresponding definitions. (Cf. Equations (36) and (37).) 

Equations (43) and (63) can be used together to develop a state-space equation for x b . 
First, express Equation (63) in measure-number form: 


( %f.= [=£»■ (Eq- 66 ) 

Next, define rotation matrix S,F Q by 

r > 

l 

= S,F Q\L 2 h (Eq 67) 

L 

v. — J y 

where the prefix indicates the rotation of frame @ relative to frame © . Finally, observe that, for 



arbitrary vectors 
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Ci = *i/, +>'i/ 2 +*i/ 3 


r -7 = x il, +y2f 2 + z 2 / 3 » 


the cross product can be expressed in determinant form by 

Iff 

£_1 £-2 i-3 

Ci x r 2 = x, .y, z, , 

*2 ^2 *2 

or in matrix form (i.e., using measure numbers) [8], by 

0 ~ z \ yi 1 1*2 

(F) ( \ 

V— i x L 2 ) = 0 -x, j y 2 

-y\ *. o 1*2 


Represent the above skew-symmetric matrix by {F) r * . Using this notation. Equation (43) can be 
expressed as follows: 

^ = (Eq 72) 

where I is the inertia matrix corresponding to 7 . 

Linearizing Equation (3) about % F = 0 yields 


N cc F = r l -M. 


(Eq. 73) 


v = 2* FI % + «, ri %. 


(Eq. 74) 


If ^ U.$ is assumed to vary negligibly with time and Equation (74) is linearized about <j> = 0 , the 


result is 


v= v% 


(Eq. 75) 


From Equations (59), (73), and (75), 
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m 


2'”p = r , M, 


or, equivalently, 


N,2 




= r'-M. 


(Eq. 76) 
(Eq. 77) 


Application of Equation (55) leads directly to 

4;( ,ls V- rl t)=\ I'M. 


(Eq. 78) 


In measure-number form, 


*.=\ s,F Qr' (F) m, 


or, equivalently, x e =- S/F Q / '[ (F) M ul + (n M ur + (F) M C ] + ~ S/F Q (F h d . 


(Eq. 79) 


(Eq. 80) 


Six state equations of the system are given by Equations (17) and (37), iterating on /; six 
more, by Equations (72) and (80). The latter six are written in terms of the various forces and 
moments acting on the system, which loads have been defined in vector form by Equations (44) 
through (48). These loads can be rewritten in measure-number form and substituted into Equations 
(72) and (80), as follows. Beginning with Equation (44), the i ,h control force can be expressed as 

(S) t c =[L, , X 5/ ^V F) 4]/, (Eq. 81) 

The resultant control force becomes 


(S £ C = E (S) £ =F c u, (Eq. 82) 

1=1 

where F' c , F c , u t , and u are defined as indicated. 

Next, using Equation (64) with (45), the translational force the umbilical exerts on the flotor 
can be expressed by 

~ ~ K t±a ~C,x b = F uta x a +F ulb x b , 


(Eq. 83) 
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where 


K< = 


K) 0 


0 

0 


0 

Kf 0 
0 K] 


(Eq. 84) 


C t = 


C) 0 


0 

0 


0 

C] 0 

0 c 3 


(Eq. 85) 


and F ut „ and F utb are defined as indicated. 


The i ,h control force F‘ c exerts on the flotor a moment M' c , defined by Equation (46). 


Using the notation introduced with Equation (71), this moment can be expressed by 


(F) M‘ = 


L, B, SIF Q ln B\l, = M‘u,. 


(Eq. 86) 


The resultant moment becomes 


(F) M C = Z (F) M = KU, (Eq. 87) 

1=1 

where Kf c and A/, are defined as indicated. 

The umbilical force F ut exerts on the flotor a moment M ul , given by Equation (47). 
Substituting from Equation (83), this moment can be expressed by 

<F) M„, = (F) Lf X S/F Q T [F uta x a + F ulb x b ] ; (Eq. 88) 

or, alternatively, (F) M, ( = M uta x a + M utb x b , (Eq. 89) 

for M uta and M utb appropriately defined. 

Finally, Equation (48) expresses the moment M„,that the umbilical applies to the flotor due 
to umbilical rotational stiffness. The following equations express in measure-number form: 


(n M,r = -r IF Q T K,x d -rl'Q'C r x. = 


(Eq. 90) 
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where 


K = 


K 

0 

0 


0 0 
K) 0 , 
0 K) 



C = 


0 

0 


0 

C, 2 

0 


0 

0 , 


and M urd and M urt are appropriately defined. 

Substituting from Equations (81) through (92), Equations (72) and (80) become, 
respectively. 


+ {j n F ^~ SIF Q ^-V F J A M^x b 

+{- s!F Q cn rr F r'M^x d +(-*"£ ^ F r l M m )x e + {^F C - SIF Q ^ f J A m]u 


- (S) a + (S) a d - s/F Q (J V ^ 

— in — a — , 


F"FL 


a A 


and 


i. =\‘ IF Qr'\M,«*'+M„ A + M urd x d + M„.x,\ 


+\ SIF Qr'M'*+\ SIF Q r ,F ’& d . 


For completeness, Equation (37) can be rewritten as x d = x e . 


To include (S) r_ Ui>E a.s states, define x c by 


(0 h (S) t N ^=i c +0) h x c , 


for some high value of circular frequency co h . Taking the Laplace Transform, 




Xc~ lS) t NrF for (o«co h . 


(Eq. 91) 
(Eq. 92) 


(Eq. 93) 

(Eq. 94) 
(Eq. 95) 

(Eq. 96) 
(Eq. 97) 


so that 


(Eq. 98) 
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Now using Equations (82), (83), 87), (89), (90), and (96) with (49), 

-“A 1 " Q m c‘ r ,r'M„)i.*»{~F r ‘i'Q tn cV,i-'M^u 


(Eq. 99) 




s / r n (f) « (/)„ 


A state-space representation of the system is given by Equations (65), (93), (94), (95), and 
(99), for state vector 


—6 




*0 


(Eq. 100) 


For the small rotation angles associated with the MIM, S/F Q is approximately equal to the 3 x 3 
identity matrix, in which case the state equations have constant coefficients. Specifically, 

L=x b , (Eq. 101) 


m J 


i^ 

m 


-{ {F Y F - F J~ X M urd )x d (Eq 102) 


- {S) a + (S) a - {F V {F) a 

—in' -d —F“F U ”d 


& =a)h im F ^~ (F ^*fE r 'K^Xa +co i~„ F «»- -co h x c 

-v h ( ^'Mjx d -co h { ^ A M m )x e +0 >l±F- 


(Eq. 103) 


+a h ( %-a} h ( % E (F) g d . 


Xd = *e 


(Eq. 104) 
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*« =\ r '[ M uu>x a + + K r(l x d + M un x e ] 


+\r'M c u+-( F h d . 

2 c - 2 - rf 


(Eq. 105) 


Concluding Remarks 

This paper has presented the derivation of algebraic, state-space equations for the Canadian 
Space Agency’s Microgravity Vibration Isolation Mount. The states employed include payload 
relative translational position and velocity, payload relative rotation and rotation rate, and payload 
translational acceleration; the relative translational position and velocity states are taken across the 
umbilical. An umbilical elongation causes a restoring force due to umbilical stiffness, and an 
umbilical elongation rate causes a restoring force due to umbilical damping. Consequently, relative 
position feedback corresponds directly to a change in effective umbilical translational stiffness; and 
relative velocity feedback corresponds, similarly, to a change in effective umbilical translational 
damping. Likewise, relative rotation and rotation-rate feedback correspond to changes, 
respectively, in effective umbilical rotational stiffness and damping. Feedback of payload 
translational acceleration causes a change in effective payload mass. Thus, a cost functional which 
penalizes the chosen states produces an intuitive effect on system effective stiffness, damping, and 
inertia values. 

The acceleration states can be selected to pertain to any arbitrary point on the flotor. This 
allows an optimal controller to be developed which penalizes directly the accelerations of any 
significant point of interest, such as the location of a crystal in a crystal-growth experiment. 

The equations have been put into state-form so that the powerful controller-design methods 
of optimal control theory (e.g., H 2 synthesis, H m synthesis, mu synthesis, mixed-mu synthesis, and 
mu analysis) can be used. The controller design approach detailed in references [9], [10], and [11] 
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has been successfully adapted for MIM controller design; the results will be presented in subsequent 
papers. 
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